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ABSTRACT 

We present the results of ultraviolet echelle spectroscopy of a sample of fifty nine 
F,G,K and M stars from the Anglo- Australian Planet Search target list. Call activity 
indices, which are essential in the interpretation of planet detection claims, have been 
determined for these stars and placed on the Mount Wilson R-hk system. 
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1 INTRODUCTION 

Extensive planet search programmes are now operating on 
a variety of telescopes around the world, including the 3.9m 
Anglo- Australian Telescope (Tinney et al. 2001, Butler et al. 
2001), the 10m Keck I (Vogt et al. 2000), the Lick 3m (Fis- 
cher et al. 2001), the Leonard Euler 1.2m (Naef 2001) and 
ESO Coude Auxiliary Telescopes (Kiirster 2000) at La Silla, 
and the Observatoire de Haute-Provence 1.9m (Queloz et al. 
1998). These programmes reach velocity precisions ranging 
from 15 to 3 ms" 1 . At these levels apparent velocity vari- 
ability (known as jitter), which is intrinsic to the target star 
itself (rather than the presence of an unseen companion), 
can be produced by the effects of stellar activity and rapid 
rotation combined with surface inhomogeneity (Saar et al. 
1998, Saar & Fischer 2000, Santos et al. 2000). 

The level of stellar activity in solar-type stars is com- 
monly examined by measuring the strength of chromo- 
spheric emission features in the cores of the Call H and 
K absorption lines. This emission is usually parametrised 
by the Mount Wilson HK Project's Rjj K index (see Bali- 
unas et al. (1998) and references therein). This index has 
been shown to be a useful indicator of the levels of jitter 
in F-M type stars (Saar, Butler & Marcy 1998). Estimates 

* Based on observations obtained at the Anglo— Australian Tele- 
scope, Siding Spring, Australia. 



of Rhk are therefore essential to planet search programmes, 
since they indicate the level of velocity variability in a given 
star which could be erroneously ascribed to the presence of 
an unseen planet, rather than jitter. 

Unfortunately, while extensive Rjjk data are available 
in the Northern Hemisphere from the Mount Wilson HK 
Project programme (Baliunas et al. 1998) and the Vienna- 
KPNO Call H&K Survey (Strassmeier et al. 2000), the only 
large Rjjk data set available in the South is that of Henry et 
al. (1996). Henry et al. measured chromospheric Call emis- 
sion in more than 800 southern stars (mostly G dwarfs). 
However, a large fraction of the stars targeted by the Anglo- 
Australian Planet Search (Tinney et al. 2001, Butler et al. 
2001) were not observed by Henry et al. We have there- 
fore begun a programme of ultraviolet (UV) spectroscopy of 
these stars on the AAT to fill this gap. 



2 OBSERVATIONS 

The Anglo-Australian Planet Search (AAPS) is being car- 
ried out on the 3.9m Anglo- Australian Telescope (AAT), 
using the University College London Echelle Spectrograph 
(UCLES) . The observations reported here were obtained on 
2001 August 4. For these observations UCLES was operated 
in its 31 lines mm -1 mode, with its UV-optimised collimator, 
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Figure 1. Example Call HK spectra from UCLES. The left panel shows an observation of HD1835, indicating the locations of the 
Call H&K lines, and the continuum V and R regions, described in Section p,2| . The fluctuations seen in this spectrum are due to weak 
absorption lines, not noise. The right panel shows an expanded region of the Call H line for an activity sequence from HD1835, HD16673, 
HD17051 to HD23429 (from top to bottom - see Table ph. Also shown are the Mount Wilson triangular Call bandpass we have usedfsolid 
line), and the Henry et al. (1996) broad Call bandpass (dotted line). All spectra have been normalised to F„=1.0 in the R, continuum 
region. The spectra have been offset by 0.4,0.1,0.2,0.0 (from top to bottom) for clarity. 



to provide an echellogram centered at 390 nm. The detector 
used was the AAO's EEV 2048x4096 13.5/xm pixel detector 
(rather than the MITLL2 2048x4096 15^m pixel CCD used 
for all previous A APS observations). This CCD provides ex- 
cellent UV response, with a quantum efficiency at the Ca II 
HK lines of ~ 65 per cent. Compared to the MITLL2 the 
EEV has smaller pixels and reduced charge diffusion, re- 
sulting in greatly improved resolution - it will become the 
default detector for all future AAPS observing. CCD gain 
was 2.7 e~ adu -1 , and read noise was 4.4 e~ per pixel. The 
CCD was binned by two spatially, to provide an effective 
slit length of 11 pixels. In the configuration adopted the dis- 
persion at the Call Rhk lines is 0.002 nm pixel -1 , while the 
resolution is 3.5 pixels or 0.007 nm. 

Sixty seven F, G, K and M-type stars were observed 
in Of! 7 seeing though a l'.'0x3'/5 slit, with exposure times 
ranging from 30s to 180s, providing a signal-to-noise ratio 
per 0.002 nm wavelength pixel of between 10 and 50. The 
stars observed included AAPS target stars which were also 
observed by Henry et al. (1996), and stars from the Mount 
Wilson HK Project, which Henry et al. used as standards for 
calibrating their observations onto the Mount Wilson Rhk 
system. The stars observed are listed in Table ^ along with 
magnitude, colour and spectral type data for each star from 
the HIPPARCOS catalogue (ESA 1997). 



3 ANALYSIS 



3.1 Image Processing 

Reduction of the CCD echellogram images proceeded in a 
standard manner in the Figaro environment. Images were 
overscan subtracted, trimmed and then flat-fielded using 
flat-field exposures from which the spectrograph blaze func- 
tion for each order had been fitted and divided out. Scat- 
tered light in the spectrograph produces a faint overall level 
of illumination of the detector - this was fitted in the gaps 
between orders and subtracted over the whole detector. 
These images were then straightened to make the echelle or- 
ders align with columns on the detector. The entire 3'.' 5 slit 
length was then extracted from each order - no sky subtrac- 
tion was performed. In the wavelength range of interest the 
sky is featureless. An observation of blank sky obtained well 
after astronomical twilight reveals sky counts contribute less 
than 0.016 e~s _1 per 0.002 nm wavelength pixel, which is 
negligible compared to our targets («8 e~ s~ per 0.002 nm 
wavelength pixel at V=7.5). 

The orders containing the Ca II H and K lines, as well as 
the adjoining orders (orders 146-142), were then extracted 
for individual processing. All were wavelength calibrated us- 
ing a ThAr spectrum obtained at the end of the night. All 
these data were also flux calibrated using an observation 
of the HST flux standard /i Col (Turnshek et al. 1990) - 
note that we use this observation only to correct the blaze 
and inter-order sensitivity of our data. The night was not 
spectro-photometric. 

The spectra in order 142 (A = 398.0 - 402.6 nm) were 
cross-correlated against the observation of HD142, and the 
barycentric correction for HD142 was applied. The result is 
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Table 1. AAPS Call R' Targets. 



Calibrators 



New Data 



Star 


v 


B-V 


Sn Tvn 
°P- A yp- 


aa 
AAT 


S 6 

Henry 


°MW 


lnp-R' 

l°g %K,AAT 


lnp-R' 


HD 1835 


6.39 


0.66 


G3V 


0.319±0.001 


0.343 


0.347 


-4.49 


-4.44 


HD 3795 


6.14 


0.72 


G3/G5V 


0.153±0.001 


0.146 


0.156 


-5.06 


-5.11 


HD 7570 


4.97 


0.57 


F8V 


0.161±0.001 


0.161 


- 


-4.91 


-4.95 


HD 10180 


7.33 


0.63 


CIV 


0.161±0.002 


0.153 


- 


-4.98 


-5.04 


HD 10700 


3.49 


0.73 


G8V 


0.163±0.000 


0.172 


0.171 


-5.00 


-4.96 


HD 11131 


6.72 


0.65 


GO 


0.311+0.001 


0.310 


0.336 


-4.50 


-4.47 


HD 13445 


6.12 


0.81 


KOV 


0.297±0.001 


0.251 


- 


-4.65 


-4.74 


HD 16673 


5.79 


0.52 


F6V 


0.209±0.001 


0.225 


0.215 


-4.69 


-4.62 


HD 17051 


5.40 


0.56 


GOV 


0.240+0.001 


0.225 


- 


-4.61 


-4.65 


HD 18907 


5.88 


0.79 


G8/K0V 


0.157±0.001 


0.109 


- 


-5.05 


-5.40 


HD 20201 


7.27 


0.58 


GOV 


0.174±0.001 


0.164 


- 


-4.88 


-4.93 


HD 20766 


5.53 


0.64 


G2V 


0.262±0.001 


0.245 


- 


-4.60 


-4.65 


HD 20782 


7.36 


0.63 


G3V 


0.180+0.001 


0.176 


- 


-4.88 


-4.91 


HD 20807 


5.24 


0.60 


G1V 


0.176±0.001 


0.195 




-4.88 


-4.76 


HD 23079 


7.12 


0.58 


F8/C0V 


1 55+0 001 


0.164 




-4.99 


-4.94 


HD 23249 


3.52 


0.92 


K0IV 


141+0 000 

yj • J. J — 1_ yj . yjyj yj 


0.129 


0.137 


-5.17 


-5.22 


HD 25874 


6.74 


0.67 


G5IV-V 


164+0 001 

\j • J. yj . yjyj _i_ 


0.173 




-4.96 


-4.93 


HD 38283 


6.69 


0.58 


C0/G1 V 


1 53+0 001 

yj • J. ^ j 'j yj i vui 


0.158 




-5.01 


-4.97 


HD 38973 

llU Ouc/ I ■ ) 


6.63 


0.59 


G2V 


n 1 >ss+0 002 

y> • ± <j o _l yj . uui 


0.175 




-5.01 


-4.88 


HD 202628 


6.75 


0.64 


G5V 


0.250+0.001 


0.217 




-4.63 


-4.73 


HD 205536 


7.07 


0.76 


G8V 


n i fio+0 001 

yj. J_»jii_l_L/. UUi 


0.179 




-5.01 


-4.94 


HD 212168 


6.12 


0.60 


G3IV 


1 57+0 001 


0.173 




-4.99 


-4.89 


HD 222335 


7.18 


0.80 


K1V 


233+0 003 


0.256 




-4.78 


-4.73 


HD 142 


5.70 


0.52 


GUV 


0.155+0.001 


- 


- 


-4.95 


- 


HD 2039 


9.00 


0.66 


G4V 


0.176+0.005 


- 


- 


-4.91 


- 


HD 2587 


8.46 


0.75 


G6V 


0.151+0.003 


- 


- 


-5.08 


- 


HD 3823 


5.89 


0.56 


G1V 


0.155+0.001 


- 


- 


-4.98 


- 


HD 6735 


7.01 


0.57 


F8V 


0.166+0.002 


- 


- 


-4.91 


- 


HD 7199 


8.06 


0.85 


K0IV/V 


0.180+0.004 


- 


- 


-4.98 


- 


HD 9280 


8.03 


0.76 


G5 


0.150+0.003 


- 


- 


-5.08 


- 


HD 10647 


5.52 


0.55 


F8V 


0.204+0.001 


- 


- 


-4.72 


- 


HD 11112 


7.13 


0.64 


G4V 


0.156+0.001 


- 


- 


-5.02 


- 


HD 16417 


5.78 


0.65 


G1V 


0.147+0.001 


- 


- 


-5.09 


- 


HD 19632 


7.29 


0.68 


G3/G5V 


0.362+0.002 


- 


- 


-4.42 


- 


HD 20029 


7.05 


0.56 


F7V 


0.149+0.001 


- 


- 


-5.02 


- 


HD 22104 


8.32 


0.68 


G3V 


0.154+0.004 


- 


- 


-5.04 


- 


HD 23127 


8.58 


0.69 


G2V 


0.163+0.004 


- 


- 


-4.99 


- 


HD 23484 


6.99 


0.87 


K1V 


0.510+0.002 


- 


- 


-4.43 


- 


HD 24112 


7.24 


0.56 


F8V 


0.152+0.001 


- 


- 


-5.00 


- 


HD 25587 


7.40 


0.54 


F7V 


0.158+0.002 


- 


- 


-4.95 


- 


HD 26754 


7.16 


0.55 


F7/F8V 


0.152+0.001 


- 


- 


-5.00 


- 


HD 27442 


4.44 


1.08 


K2IV 


0.150+0.001 


- 


- 


-5.27 


- 


HD 28255 


6.28 


0.66 


G4V 


0.181+0.001 


- 


- 


-4.89 


- 


GJ9155B d 


7.50 


0.71 


G6V 


0.258+0.001 


- 


- 


-4.65 


- 


HD 30177 


8.41 


0.77 


G8V 


0.151+0.004 


- 


- 


-5.08 


- 


HD 30876 


7.49 


0.90 


K2V 


0.469+0.002 


- 


- 


-4.51 


- 


HD 31527 


7.49 


0.61 


G2V 


0.165+0.001 


- 


- 


-4.95 


- 


HD 36108 


6.78 


0.59 


G3V 


0.153+0.002 


- 


- 


-5.01 


- 


HD 38382 


6.34 


0.58 


F8/G0V 


0.162+0.001 






-4.94 




HD 40307 


7.17 


0.94 


K3V 


0.267+0.004 






-4.83 




HD 204385 


7.14 


0.56 


G0IV 


0.158+0.001 






-4.98 




HD 204961 


8.66 


1.52 


M1V 


0.975+0.008 






-5.10 




HD 205390 


7.14 


0.88 


K2V 


0.424+0.002 






-4.53 




HD 209268 


6.88 


0.56 


F7V 


0.153+0.001 






-4.99 




HD 211317 


7.26 


0.65 


G4IV 


0.156+0.002 






-5.03 




HD 216437 


6.04 


0.66 


G4IV-V 


0.154+0.001 






-5.04 




HD 217987 


7.34 


1.48 


M2/M3V 


1.074+0.004 






-5.01 




HD 222237 


7.09 


0.99 


K3V „ 


0.281+0.007 






-4.86 





lno- R ' c 

-4.44 
-5.04 



-4.96 
-4.43 



-4.66 



-5.18 



Notes : a - Saat calculated as described in Section g after calibration onto the Mount Wilson Smw index, b - Suenry is the Call HK 
index of Henry et al. (1996) (calibrated onto the same scale as the Mount Wilson Smw index), as provided in columns 8-11 of their 
Table 3, and column 14 of their Table 4. Smw is the Mount Wilson index used to perform this calibration from columns 4-7 of their 



Table 3. 



5R! 



HK.CTIO 



are the logRC, K indices measured by Henry et al. (1996), while logRj, 



are the log R!r K indices Henry et 



al. used to check their calibration, d - Also known as SA0258997. This star lies within 6" of HD28255 (itself also known as GJ9155A) 
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a set of zero- velocity spectra. (Though HD142 is a radial 
velocity variable at the 30m s _1 level (Tinney et al. 2002), 
this variability is negligible for our purposes). Lastly, all the 
spectra were re-sampled to the same 0.002 nm binning. 



3.2 S Indices 

The resulting spectra were then used to produce a chromo- 
spheric Call HK emission index, Smw- The Mount Wilson 
HK Project defines its Smw (Duncan et al. 1991) as the ra- 
tio of the flux in two triangular bands with full-width at 
half-maximum of 0.109 nm centered at 393.3664 (K) and 
396.8470 nm (H), to that in two 2.0 nm wide continuum 
bands centered at 390.1 (V) and 400.1 nm (R). This ratio 
is measured using a specifically-designed multi-channel pho- 
tometer. Smw defined as 

N H + N K 



Nr + Nv' 



where the Ni are the counts in each of the H, K, V and R 
bands, and a is a constant determined for each night by the 
observation of standards. 

Using the high-resolution spectra obtained at the AAT, 
we were able to define our Call H and K and continuum 
bands as having the same effective profiles (i.e. 0.1 nm tri- 
angular and 2.0 nm square bandpasses) as used at Mount 
Wilson. At the AAT, therefore, we measure the flux in each 
band, and normalise for the band width to obtain 



N : 



E A f(A)AA 



in each band. Because both the Ca H and Ca K lines appear 
in two orders in our echellogram, we actually obtain two flux 
estimates for each, which we combine to form an observed 
Saat index as follows 

„ = N H ,i + N H ,2 + Aam + N K:2 
AAT ~ 2{N R + N V ) 

In Figure ^| we plot this observed Saat against the Smw 
and SHenry values reported by Henry et al. A tight linear 
correlation is seen between Smw and Saat, with rms scat- 
ter about a linear fit of 0.02, zero-point offset of 0.061 and 
slope of 1.043. Examination of Fig. ^ would suggest that 
a calibration against the MW data alone would have even 
smaller scatter, and could be preferred. Certainly our high- 
resolution spectra should transfer onto the MW system in 
a linear fashion with no systematic effects. Unfortunately, 
the number of MW stars accessible in the south is small, so 
such a calibration would (at present) be based on a small 
number of stars, and only a few active ones. Because such 
active stars are by their very nature highly variable in their 
activity, we have therefore chosen to use all the MW and 
Henry et al. calibrators to derive the relation shown in Fig. 
^| No major correlation is seen between degree of deviation 
from the fit and spectral type. 

The values of Saat so calibrated (including propagated 
photon-counting uncertainties), and the calibrating S in- 
dices are shown in Table |l[ In general the scatter about 
the calibration curve is significantly larger than the photon- 
counting uncertainties, which likely reflects the intrinsically 
variable nature of stellar activity. 

It should be noted that an advantage of performing 




0.4 



Figure 2. Calibration of the Saat spectral index on to the S 
index scale of Mount Wilson (open symbols) and Henry et al. 
(1996) (filled symbols) surveys, for the observed calibration ob- 
jects. Different spectral types are plotted with different symbols. 
The dotted line shows a unitary transformation, and the solid line 
is the fitted linear calibration (Smw = 1-043 Saat + 0.061). 



these observations at high spectral resolution is that we can 
define our bandpasses in exactly the same way as those used 
at Mount Wilson - the resulting indices require only a linear 
calibration with slope of 1.0 to place them on a standard sys- 
tem. Henry et al., in contrast, were forced to use a non-linear 
calibration, which changes from run to run, introducing pos- 
sible systematic effects. Experiments on our data show that 
when the Call bandpass is broadened, non-linear calibra- 
tions of increasing complexity become necessary, with the 
consequent risk of introducing systematic uncertainties. 

3.3 R-hk Indices 

The S index provides an estimator of photospheric plus 
chromospheric Call strength. For activity purposes, how- 
ever, it is the chromospheric emission which is of interest. 
This is usually normalised to the star's bolometric lumi- 
nosity, and parametrised by the Rhk index. The derivation 
of Rhk from S is described in detail in the Appendix of 
Noyes et al. (1984). We follow these procedures to derive the 
log Rhk. aat values shown in Table hi where we also show: 
the logR HK CT i values measured by Henry et al. (1996); 
the log R HK mw values used to calibrate those data; and the 
HIPPARCOS B-V photometry used to estimate the photo- 
metric contribution to Saat. 

In Figure || we plot log Rhk, aat against log Rhk,mw 
and logR HKCT i - Note that there are no arbitrary cali- 
bration factors included in the transformation from Saat 
to Rhk aat- The resulting plot is consistent with a uni- 
tary transformation - the solid line fit shown in Fig. |^ is 
for a zero-offset of 0.09±2.3 and a slope 1.02±0.09, and 
has rms scatter of 0.09 (for all data), or 0.05 for the MW 
data alone . We therefore apply no further calibration to 
our log R H k, a at estimates, presented in Table but note a 
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Figure 3. Plot of the derived R H k AAT versus the Mount Wil- 
son (open symbols) and Henry et al. (1996) (filled symbols) R H k 
measurements for our calibration objects. The dotted line shows 
a unitary transformation, and the solid line is a linear fit. 

root-mean-square scatter about the MW standards of ±0.05 
(in the logarithm), which is a better uncertainty estimate 
for the final results than the photon-counting uncertainties. 
The rms scatter about the fit for observations in common 
with Henry et al. is larger at ±0.09. Expressed as the "mean 
absolute deviation" (used by Henry et al.) this corresponds 
to ±0.07, while the mean absolute deviation for the calibra- 
tion of the Henry et al. sample onto the MW data is ±0.052, 
leading us to conclude a large part of this scatter may be 
intrinsic to the Henry et al. comparison sample rather than 
the AAT data. 



4 DISCUSSION 

Table |l| provides at least two new Rhk measurements for 
stars with known planetary companions, HD27442 (Butler et 
al. 2001) and HD142 (Tinney et al. 2002) - though there may 
be several more, as yet undiscovered, planets in this set of 
stars (for which AAPS observing is on-going). In the case of 
HD27442, we find the star, at this epoch, to be very inactive. 
Saar et al. (1998) and Santos et al. (2000) have examined 
the levels of velocity variability intrinsic to stars themselves 
(a' v ), as a function of activity (Rhk)> m stars monitored 
by the Lick 3m planet search. A trend emerges in which 
more active stars show more intrinsic velocity instability, or 
"jitter", with a form of a' v oc(Rhk) 1 ' 1 for G and K dwarfs, 
and with a' v in the range 3-10 ms~ for log Rhk =-5.0. 

Based on this we would predict the jitter in HD27442 
to be £T.5-5 ms -1 . Butler et al. (2001) derive a scatter 
about their Keplerian fit to the velocity data for HD27442 of 
2.9 ms -1 , which is an upper limit to the level of jitter which 
could be present in this star - in line with its extremely low 
level of Call activity. 



HD142 is only slightly more active at logR HK =-4.95, 
implying an expected level of jitter in the range ^4-12 m s _1 . 
Tinney et al. (2002) report a scatter about a single Keplerian 
fit of 5.9 ms -1 , consistent with this prediction. 

Several of the stars in Table ^ reveal themselves to be 
quite active, with logR HK > -4.5. This would correspond to 
likely levels of jitter greater than 10-30 m/s in HD19632, 
HD30876 and HD205390. Planet candidates detected around 
stars like this must be scrutinised with great care to avoid 
erroneous detection claims. The recent results of Queloz et 
al. (2001) for HD166435, in particular, provide a valuable 
cautionary tale. HD166435 was selected for observation as 
part of the Observatoire de Haute- Provence ELODIE planet 
search. The available data (which did not include Rhk mea- 
surements) did not indicate particular youth or activity. 
Over the period May 1997 to September 1999 seventy spec- 
tra of this star were acquired, as it appeared to have an 
unseen planetary companion. Further analysis, however, re- 
vealed this variability to be entirely due to activity. A result 
which is not unexpected when the value of the S index of 
this star (as measured at Mount Wilson) is considered. At 
S=0.42-0.51 (Queloz et al. 2002) this GO dwarf is phenome- 
nally active - more active than the most active stars shown 
in Figs. ^ and [| implying logR HK ~ —4. At this level of 
activity the Saar et al. (1998) results would imply an ex- 
pected level of jitter in the range 36-120 ms -1 , which would 
go much of the way toward explaining the 200 ms -1 am- 
plitude velocity variations seen in HD164435. Having R H k 
measurements in hand when considering the validity of a 
planet candidate is clearly essential. 

Finally we would note the intrinsically variable nature 
of active stars, and the importance of obtaining Ca II obser- 
vations at multiple epochs in order to be sure a star really 
is inactive. The reliability of such data is important when 
trying to decide what level of velocity variability could be 
due to activity /rotation induced jitter, and what level can 
be trusted to be due to an unseen companion. 



5 CONCLUSION 

Call HK activity measurements are provided for 59 stars 
from the Anglo- Australian Planet Search target list - thirty 
five of which are for stars not previously published. Four- 
metre telescope UV echelle spectroscopy is an efficient way 
in which to obtain high quality data for this purpose, and 
because of the high resolution of such data, calibration onto 
the standard Mount Wilson R H k system can be carried out 
with greatly reduced systematic uncertainties. 
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